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Abstract 

We start with some global Maxwellian function M, which is a sta- 
tionary solution (with the constant total density p) of the Boltzmann 
equation, and we denote the number of the corresponding space vari- 
ables by n. The notion of distance between the global Maxwellian 
function and an arbitrary solution / (with the same total density p at 
the fixed moment t) of the Boltzmann equation is introduced. In this 
way we essentially generalize the important Kullback-Leibler distance, 
which was used before. An extremal problem to find a solution of the 
Boltzmann equation, such that dist{M, /} is minimal in the class of 
solutions with the fixed values of energy and of n moments, is solved. 

1 Introduction. 

The well-known Boltzmann equation for the monoatomic gas has the form 

^ = -C-v,/ + Q(/,/), (1.1) 
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where ieR stands for time, x = (xi, ...,x n )EQ stands for space coordinates, 
( = (d, C n )GM n is velocity, and K denotes the real axis. The collision 
operator Q is defined by the relation 

Q(fJ)= I I [fiOfiQ-fiOfiQWC-c^dadC, (i.2) 

where B(( — cr)>0 is the collision kernel. Here we used the notation 

C = (C* + 0/2 + or|C* - CI/2, C = (C* + 0/2 - a|C - CI/2, (1.3) 

where <7GS' ri_1 , that is, erGM™ and \a\ = 1. The solution f(t,x,Q of Boltz- 
mann equation (11.11) is the distribution function of gas. We start with some 
global Maxwellian function M, which is the stationary solution (with the 
total density p) of the Boltzmann equation. The notion of distance between 
the global Maxwellian function and an arbitrary solution / (with the same 
value p of the total density at the fixed moment t) of the Boltzmann equa- 
tion is introduced. In this way we essentially generalize the Kullback-Leibler 
distance [J], which was fruitfully used before (see further references in the 
recent papers [2],[9j[T2]). Our approach enables us to treat also the non- 
homogeneous case. An extremal problem to find a solution of the Boltzmann 
equation, such that dist{M, /} is minimal in the class of solutions with the 
fixed values of energy and of n moments, is solved. 

Some necessary preliminary definitions and results are given in Section 
An important functional, which attains maximum at the global Maxwellian 
function is introduced in Section [3j The distance between solutions and the 
corresponding extremal problem are studied in Section |H 

We use the notation Cq to denote the class of different iable functions 
/(C), which tend to zero sufficiently rapidly when ( tends to infinity. 

2 Preliminaries: main definitions and results 

In this section we present some well-known notions and results connected 
with the Boltzmann equation. The distribution function f(t,x,() is non- 
negative: 

f{t,x,C)>0, (2.1) 

and so the entropy 

S(t,f) = -[ f f(t,x,()logf(t,x,0d(dx (2.2) 
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is well-defined. 

Definition 2.1 A function <p(Q is called a collision invariant if it satisfies 
the relation 

[ 0(C)W,/)(CR = O for all feC 1 ,. (2.3) 

JM. n 

It is well-known (see [TT]) that there are the following collision invariants: 
<MC) = 1, MC)=Q (i = l,2,...,n), n+1 (C) = |C| 2 . (2.4) 

The notions of density p(t, x), total density p(t), mean velocity u(t, x), energy 
E(t,x), and total energy E(t) are introduced via formulas: 



p(t,x) = J f(t,x,()d(, p(t) = j p(t,x)dx, (2.5) 
u(t,x) = {l/p(x,t)) J tf(t,x,()d(, (2.6) 

E(t,x)= /f/(M,CK, E(t)= [ I l -^f(t,x,()d(dx. (2.7) 
The function 

/ = (p/(2vrT)"/ 2 ) exp ( - |C - u\ 2 /(2T)). (2.8) 

is called the global Maxwellian and is a function of the mass density p > 0, 
bulk velocity u = (u±, ...,u n ) and temperature T. We assume that the domain 
Q is bounded and so its volume is bounded too: 

Vol(fi) = V n < oo. (2.9) 

Therefore, the function 

M(C) = (p/(V^(2vrT)"/ 2 )) exp ( - |( - u\ 2 /{2T)) (2.10) 
is a global Maxwellian with the constant total density p. 

Proposition 2.2 [77]/ ITte global Maxwellian function M(() is the station- 
ary solution of the Boltzmann equation (11. ip . 

Boltzmann proved in pQ the fundamental result below: 

Theorem 2.3 Let f G Cq be a non-negative solution of equation ( II. ip . T/ien 
i/ie following inequality holds: 

dS/dt > 0. (2.11) 
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3 Extremal problem 



Similar to the cases considered in [7J[8], an important role is played by the 
functional 

F(f) = XE + S, A = -l/T, (3.1) 

where S and E, respectively, are defined by formulas (12.21) and (12 .7p . The 
parameters A = —1/T and p are fixed. Now, we use the calculus of variations 
(see [3]) and find the function / max which maximizes the functional (13. ip on 
the class of functions with the same p(t) = p at the fixed moment t. The 
corresponding Euler's equation takes the form 



5_ 
Sf 
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\\^-f-f\ogf + pf =0. (3.2) 



Here jj stands for the functional derivative. Our extremal problem is con- 
ditional and /i is the Lagrange multiplier. Hence, we have 

A ^!_l_l og/ + /i = 0. (3.3) 

From the last relation we obtain 

/ = Ce -|ClV(2T)_ (34) 

Formulas (I2.10p and (13.41) imply that 

f = M(Q = ——^-—e--£. (3.5) 

y 1 Vh(27rT)"/ 2 v ' 



We have the inequality 

6 2 



_F=-l//<0. (3.6) 



Corollary 3.1 The global Maxwellian function M((), which is defined by 
formula (13. 4p . gives the maximum of the functional F on the class of func- 
tions with the same value p of the total density pit) at the fixed moment 
t. 



4 



In view of O, flUD, and (EHJ we see that 

(F(f))(t) = -[ [ (^ + \ogf(t,x,C))f(t,x,Od(dx. (3.7) 
It follows from (|2~5]1 . f l33|) . and (ETTJ) that 

F(M) = ^ l0g (M2^)- (3 ' 8 > 

Therefore, Corollary 13.11 can also be proved without using the calculus of 
variation (see [ID]). Indeed, taking into account relations (I3.5p . (13 .7p . and 
(13.81) and the fact that the total densities of M and / are equal, we have 

F(M) - F(f) = 1 1 M(l - L + L loe L) d(dx . (3.9) 

Using inequality 1 — x + a; log a; > for x > 0, Xy^l, we derive from (13.91) 
that 

F(M)-F(f)>0 (f^M). (3.10) 



Remark 3.2 Since the extremal problem is conditional, the connection be- 
tween the energy and entropy can be interpreted in terms of game theory. The 
functional (13. ip defines this game. The global Maxwellian function M(() is 
the solution of it. A game interpretation of quantum and classical mechanics 
problems is given in the papers ^\SI- 

4 Distance 

Let f(t, x, () be a nonnegative solution of the Boltzmann equation (II. ip . We 
assume that T and the value p = p(t) at some moment t are fixed. According 
to (13.101) we have 

F(M)-F(f)>0, (4.1) 

where the global Maxwellian function M(() is defined in (13. 5p . The equality 
in (14. ip holds if and only if f(t,x,() = ^(0- Hence, we can introduce the 
following definition of distance between the solution f(t, x, () and the global 
Maxwellian function M((): 

tist{M,f} = F{M)-F(f). (4.2) 
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Remark 4.1 In the spatially homogeneous case, if not only the total den- 
sities pm and pj of M and f are equal, but the energies Em and Ef are 
equal too, then our definition (14. 2 p of distance coincides with the Kullback- 
Leibler distance (see IHtf ). However, our approach enables us to treat also 
the non-homogeneous case. 

Next, we study the case Em ^ Ef and start with an example. 
Example 4.2 Let T x ^ T and consider the global Maxwellian function 

M ^ = v^wr^{ J id- < 43 > 

Direct calculation shows that 

E x = E Ml = pnT x /2 + E, (4.4) 
F(M X ) = -p(log { y^ Ti)n/2 ) ~ n(l ~ T x /T)/2). (4.5) 

It follows from (ESD and (gj2) that 

dist{M, Mi} = -pn(\og{T x /T) - T x /T + l)/2. (4.6) 

We introduce the class C(p, E x , U) of non- negative functions functions f(t, x, () 
with the given total density p (see (2.8)), total energy 

/ / ^f(t,x,()d(dx = E x , (4.7) 

and total moments U — (U x , U%, U n ) , where 



Uk= ( k f(t,x,C)d(dx. (4.8) 

Recall that the global Maxwellian function M is defined by (13. 5p . 

Extremal problem. Find a function f , which minimizes the functional 
dist{M, /} on the class C(p, E X ,U). 

The corresponding Euler's equation takes the form 



d r (A + v)^-f - /log/ + pf + /^ 7 fcCfc] =0- (4.9) 

k 
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Recall that our extremal problem is conditional, and p, v, 7^ are the La- 
grange multipliers. Hence, we have 

(A + u)^- - log / - 1 + fi + J2^Ck = 0. (4.10) 

k 

From the last relation we obtain 

ICI 



/ = Cexp ((A + + 7*Cfc) • (4-H) 

k 



According to (12. 5p we have A + v < 0. Now, we rewrite (14. lip as 

/^-^r^^E(^) 2 ). 

where 

To calculate the parameters //, z/, 7^ we use again the well-known formulas 

/DO /»00 1 

e" ^ = / € 2 e- a ?d£=—y/n/a, a > 0. (4.14) 

-OQ J— OO ^ 

Formulas (12511 . (H77I) . (Oj) . f l4TT2|) . and fl4TT4j) imply that 

C! = p/V n , lk /(X + v) = -U k /p, -(\ + u) = 7? 1 , (4.15) 

where 

np np z / — ' 

Because of (I4.12p and (14.151) we see that / is just another global Maxwellian 
function 

/ = M >«>=T«^-K- K ^r £ )- < 4 - 17 > 



77j[dist{M,/}] = l// (4.18) 



Moreover, the inequality 

i? 

holds, that is, the functional dist{M, /} attains its minimum on the function 
/ = Mi, which satisfies conditions pit) = p, ( 14. 7ft . and ( 14.8D . The following 
assertion is true. 
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Proposition 4.3 Let M and Mi, respectively, be defined by (j3.5p and (14.171) . 

If the function f satisfies conditions p(t) = p, (14. 7p . (14. 8p . and f^M\, then 

dist{M,/}>-^(log(T 1 /T)-r 1 /T + l)+ |( h 



2 v ov ' ' ' ' 2pT } 
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